Quantum state tomography of an itinerant squeezed microwave field 
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We perform state tomography of an itinerant squeezed state of tiie microwave field prepared by 
a Josephson parametric amplifier (JPA). We use a second JPA as a pre-amplifier to improve the 
quantum efhciency of the field quadrature measurement (QM) from 2 % to 36 ± 4 %. Without 
correcting for the detection inefhciency we observe a minimum quadrature variance which is 681? % 
of the variance of the vacuum. We reconstruct the state's density matrix by a maximum likelihood 
method and infer that the squeezed state has a minimum variance less than 40 % of the vacuum, 
with uncertainty mostly caused by calibration systematics. 
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Fundamental quantum optics experiments at mi- 
crowave frequencies have been recently performed with 
superconducting qubits or Rydberg atoms inside high- 
quality microwave cavities. Examples include the recon- 
struction of the Wigner functions of Fock states from one 
[H to a fewphotons and coherent superpositions of few 
photons [24J1- States such as these, which are manifestly 
nonclassical light states, are crucial for quantum infor- 
mation processing, because they can be used to generate 
entanglement. However, in the cited experiments, these 
states are confined in cavities. Therefore distributing 
entanglement to separate parties, as required in quan- 
tum communication protocols, remains challenging for 
microwave implementations. In contrast to the discrete 
Fock state approach, continuous variables quantum infor- 
mation (CVQI) strategy uses another type of nonclassi- 
cal states, the squeezed states, which are readily created 
in itinerant modes. These states exhibit reduced noise, 
below the vacuum fluctuations, in one of their quadra- 
ture components and amplified noise in the other one. 
They are also easily generated at optical frequencies in 
the itinerant output modes of parametric amplifiers made 
of optically nonlinear crystals. At optical frequencies, 
CVQI has progressed rapidly from the initial creation of 
squeezed states Q and tomographic reconstruction [^-Q 
of those states to teleportation and quantum error 
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At microwave frequencies, the field is less advanced. 
The generation of microwave squeezed states using the 
nonlinear electrical response of superconducting Joseph- 
son junctions has been reported [l3|, with inferred 
squeezing down to 10 % of vacuum variance [l^ . Such 
states can be powerful tools for quantum information 
processing and communication because microwaves and 
superconducting qubits can mimic useful light-atom in- 
teractions, as demonstrated in [l5[. Furthermore, these 
devices are made of compact and integrable electrical cir- 



cuits, with much promise for building complex quantum 
information processors. The lack of an efficient quadra- 
ture measurement (QM) for itinerant modes has slowed 
the advancement of CVQI. However, as demonstrated re- 
cently in 1 it is possible with a JPA to realize an effi- 
cient single QM. 

In this Letter, we report the tomography of an itinerant 
squeezed microwave field. We demonstrate that our JPA 
based measurement scheme has a quantum efficiency 20 
times greater than a QM employing statc-of-art semicon- 
ductor amplifiers. We infer the quantum state prepared 
by maximum likelihood tomography, correcting for inef- 
ficiency in our QM. We discuss the achieved degree of 
squeezing, from the perspective of generating entangle- 
ment on chip. 

Homodyne tomography is a standard experimental 
tool to infer the quantum state of a single mode of light. 
It was proposed in [l7| and pioneered on a squeezed op- 
tical field in Q . Its principle is depicted in the Fig. [1] A 
homodyne detection apparatus measures the value of the 
quadrature Xg, where 6 is set by adjusting the phase 
of the local oscillator. The probability density func- 
tion pT(Xg) for measuring a particular value of Xg is 
the marginal density function of the Wigner function, 
i.e. pr(A:6i) = / dXg+T,/2W{Xg, Xgj^^/2), as shown in 
Fig. [1] (b). Thus by performing measurements of Xg 
on many identical copies of the state and varying 9, the 
"hidden" quantum object can be seen from different an- 
gles and its state inferred. Losses and other Gaussian 
noise sources in the homodyne detector can be modeled 
with the insertion of a fictitious beam splitter of trans- 
missivity 77, as shown in Fig. [1] (a). In such a case, the 
measured pi{Xg) are no longer the projections of the de- 
sired Wigner function W , but of a smoother distribution 
which is the convolution of W with a Gaussian Wigner 
function [3l • However methods like maximum likelihood 
quantum state tomography can be used to deconvolve the 
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FIG. 1: Principle of the experiment, (a-left); The squeezer 
(SQ, in red) prepares a squeezed state whose quadrature dis- 
tributions are measured for different phases 6 witli an effi- 
ciency T). (a-right): Simulated measurement results for 20,000 
realizations of creating the squeezed state and measuring it 
at a single 6. The top graph shows the measured quadrature 
value versus realization number. The bottom plot is a his- 
togram (blue circles) and Gaussian probability distribution 
pxi^Xe) (red curve) of this random process, (b): Graphical 
Interpretation: the probability distribution pr(Xe) is simply 
the projection of the Wigner function. 



effect of inefficiency [l9| . 

At optical frequencies, 77 > 90 % is routinely obtained 
using a pair of balanced photodiodes [l^. Such detec- 
tors are not available for microwaves and until recently 
the best setup was a chain of phase-insensitive amplifiers 
followed by a mixer, or two such chains in parallel 
221 . In such a case, noise A„ greater than 1/2 (the vac- 
uum variance) must be added to the QM . This noise 
can be modeled as an effective efficiency by the relation 
77 = 1/(1 -I- 2 An) |2J], so the QM efficiency using phase- 
insensitive amplifiers is limited to 50 %. State of the art 
microwave amplifiers, high-electron-mobility transistors 
(HEMTs), have yl„ sa 10 — 20. In practice, the unavoid- 
able losses present in a microwave experiment t ypi cally 
result in r/ w 2%. However, as demonstrated in [l6[, in- 
serting a JPA used as a single quadrature preamplifier 
before the HEMT increases the experimentally achieved 
?7 by a factor of approximately 20. 

To perform a high-quality reconstruction of the Wigner 
function of a squeezed microwave state we operate two 
JPAs in series, as shown in Fig. [2] (b). The first JPA, 
referred to as the squeezer (SQ), prepares the squeezed 
state. The second JPA, referred to as the pre-amplifier 
(AMP), amplifies the quadrature of the squeezed state 
determined by the phase difference 6 between the AMP 
and the SQ pump tones. We vary 6 by applying to 
the two cavities pump tones slightly detuned from one 
another. The SQ stage is pumped at 7.45 GHz, while 
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FIG. 2: (a): To implement a SQ or AMP at microwaves, 
three microwaves components are required: (i) a JPA, (ii) a 
directional coupler and (iii) a circulator. Taking port (1) of 
the directional coupler as reference, (2) is the weakly coupled 
port, (3) the isolated port and (4) the direct port. Port (2) is 
used to pump the JPA. Port (3) is used to apply a cancela- 
tion tone (adjusted with a room temperature attenuator and 
phase shifter) that nulls the pump and displaces the output of 
the JPA back to the origin of the phase space, (b): Schematic 
of the experiment. In this figure, all the microwave compo- 
nents and cables are considered lossless; their imperfections 
are absorbed into the experimentally determined total trans- 
missivities ^, a and /3. 



the AMP stage is pumped at 100 kHz higher frequency; 
therefore, sweeping through 27r every 10 jis. 

Our implementation of an SQ or an AMP at mi- 
crowaves, as shown in Fig. [5] (a), requires three elements: 
(i) a JPA used in reflection, (ii) a directional coupler and 
(iii) a circulator. As described in [TJ], the JPAs are non- 
linear resonant cavities built from coplanar waveguides 
whose central conductor has been replaced by a series 
of many Josephson junctions. The Josephson junctions' 
nonlinearity causes the cavity's phase velocity to be in- 
tensity dependent. Therefore when the cavity is pumped 
it becomes a phase sensitive amplifier for input modes 
whose frequencies lie within the bandwidth of the JPA 
centered on the pump frequency. Such microwave modes 
incident on the JPA are reflected and exit the cavity 
with one quadrature amplifled and the other squeezed, 
depending on their phase relative to the pump's phase. 
A directional coupler is used to add the pump tone to 
the incident signal and remove the pump tone from the 
reflected signal. Finally the incident and reflected modes 
are separated into different cables using a circulator. 

Following Fig.[2](b), in the limit of large HEMT power 
gain Gh, our quantum efficiency can be cast as 
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2 + 2AA-a+ [2Ah - (1 - P)]IGaP ' 



(1) 



where A a {Ah) is the AMP (HEMT) added noise, a (/3) 
is the fraction of power transmitted by the microwave cir- 
cuitry between the SQ and the AMP (the AMP and the 
HEMT), and Ga is the power gain of the AMP stage. A 
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detailed description of how we calibrate each of these pa- 
rameters is in the supplementary information. Briefly, we 
inject different amounts of thermal noise into the ampli- 
fier chain while operating each JPA either as an amplifier 
(ON) or as a noiseless element with unit gain (OFF). We 
then infer the added noise and loss of the elements by 
observing the variation in the noise at the output of the 
measurement chain. The thermal noise is varied by con- 
necting the input of the SQ through a switch to either a 
"hot load" (50 microwave termination at 4.1 K) or a 
"cold load" (at 20 mK). Although the tomography is only 
performed with the "cold load" , both are required for cal- 
ibration. We obtain Aa = 0.25 ± 0.06, Ah = 17.3 ± 0.1, 
a = 68±2 % and /? = 74±5 %. However, as the switch is 
operated at the 4.1 K stage and is slightly lossy, the state 
presented at the input of the SQ with the "cold load" is 
not pure quantum vacuum, but a low occupancy thermal 
state with average photon number n ~ 0.15 ± 0.15. One 
quadrature of the resulting squeezed state is then ampli- 
fied at the AMP stage with sufficient gain Ga = 180 such 
that the noise in the amplified quadrature exceeds Ah for 
any 9. From Eq. ([T|), we obtained an overall quantum ef- 
ficiency of 36 ± 4 %, which can be compared to 77 w 2 % 
without the AMP stage. 

In this experiment our uncertainty in ?/ and n cre- 
ate a systematic source of error. We thus perform our 
data analysis under three assumptions (1) high efficiency 
(77 = 0.40) and high mean photon number {n = 0.30), 
(2) best estimate for both efficiency (77 = 0.36) and 
mean photon number (n ~ 0.15), and (3) low efficiency 
(r] = 0.33) and low mean photon number (n = 0). These 
three cases give us "pessimistic" , "best-guess" , and "op- 
timistic" analyses, in terms of the purity of the squeezed 
state estimated by the tomography. Using a lower esti- 
mate for rj and n as inputs to the tomography algorithm 
causes it to return a more pure, more squeezed, and 
therefore a more "optimistic" estimate of the squeezed 
state. Associated with each of these three cases, we also 
have statistical uncertainty, so the given error bounds 
cover an interval that includes both uncertainties around 
the "best-guess" estimate. They are reported in the form 
X^Y^ where X is the statistical mean using the "best- 
guess" calibration and L and U are respectively the lower 
and upper bounds of the one standard deviation uncer- 
tainty in the "pessimistic" and "optimistic" cases. 

We must also calibrate the QM to convert the mea- 
sured voltage noise into units of noise quanta. In optical 
homodyne tomography, this is usually done by inserting 
the vacuum and observing the quadrature noise. Anal- 
ogously, we insert the weak thermal state with mean 
photons n (by simply turning the SQ stage OFF) and 
measure voltages proportional to quadrature values at 
many 9, as shown (in blue) in Fig. |3l As expected 
this voltage noise is 9 independent, with a variance 
AVg^Q Qpp = 3.2 X 10~^ mV^. Under the convention 
that vacuum has variance 1/2 in unitless quadrature 
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FIG. 3: (a): Density plot of number of occurrences in a 1 pV 
bin size of the amplified quadrature voltage Ve versus 6/2n, 
with the SQ pump OFF (top) and ON (bottom), (b): In par- 
ticular, histograms of Ve at the maximum of squeezing: data 
(o) and Gaussian fit (continuous lines) for the SQ pump OFF 
(blue) and ON (red), (c): Noise variance AXg in quanta units 
on a log scale versus d/2n for the SQ pump ON (red) and OFF 
(blue). The (black) line indicates our estimate of the vacuum 
noise level under the "best-guess" calibration. 



space (or in units of "quanta"), we calibrate this volt- 
age variance to AXgQ qff = (1 ~ + ?7(l/2 + n) = 
0.55lo;g5 quanta. Therefore the desired conversion factor 

A^iQ,OFF/AK?Q,oFF = 1.7ltS;f7X 10^ quanta/mV-' is 
used to rescale the variances in Fig. [3] (c) . 

In Fig. [3] (a), we show QM data of the squeezed state. 
With SQ ON (red) we observe the characteristic phase 
dependent noise for a squeezed state. At the phase for 
which the variance is minimum, we show the histogram 
of quadrature measurements in Fig.[3](b). The SQ OFF 
histogram is clearly wider than the SQ ON histogram, 
demonstrating our ability to observe squeezing directly 
at the output of our measurement chain. In Fig. [3] (c) 
we plot the variance of the QM with SQ ON and OFF 
as a function of 9, expressed in units of quanta, clearly 
showing squeezing below the vacuum level. Without cor- 
recting for 7], we observe a minimum quadrature variance 

which is aa:|q j^jN 



68I7 % of the vacuum variance. 



To infer the quantum state created by the squeezer, 
correcting for loss during the QM, we used maximum 
likelihood quantum state tomography [2^. For each of 
the three calibration cases, we performed 35 reconstruc- 
tions using independent subsets each containing 10,000 
QMs of the total measured data. We estimated statis- 
tical uncertainty from the spread of properties (such as 
fidelity or minimum variance) of the set of 35 reconstruc- 
tions. The statistical uncertainty was significantly lower 
than the systematic uncertainty. In Fig. |4] we show the 
Wigner function of the "best-guess" reconstructed state 
p. The pure squeezed vacuum state \ip) that has the 
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(as both squeezers and phase-sensitive amphfiers) will re- 
quire more detailed investigation of the source of this 
noise. We used maximum likelihood quantum state to- 
mography to deconvolve the QM inefficiency in order to 
precisely characterize the state generated. This is an im- 
portant step toward generating easily distributable mi- 
crowave entanglement on chip. 

Notes: A different method was recently used to obtain 
a similar state reconstruction [2^ . 

The authors acknowledge support from the 
DARPA/MTO QuEST program. 



FIG. 4: Mean of 35 reconstructions of the Wigner function of 
the state exiting the SQ, inferred by maximum likelihood un- 
der the "best-guess" assumption, in quanta units. The faint 
pattern of ripples extending from the origin is caused by trun- 
cation at 30 photons of the density matrix used to represent 
the state. The white circle at the origin shows the full-width 
at half-maximum of the vacuum state. 



highest fidelity with p has minimum quadrature variance 
6.0^2^ % of the vacuum variance, and that maximum 



fidelity is = {ip\p\il)) 



0.81 



-1-0.16 
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As explained in 



the supplementary information, the minimum variance 
of p is biased by an amount comparable to our sys- 
tematic uncertainty, so we infer the minimum variance 
Aa;gQ i^jj^ directly from the observed minimum variance 
(l/r/)(AX|Q_MiN - (1 - We find 



as AxgQjy[jp^ 



121^2 % of the vacuum variance. For com- 



parison, the most highly squeezed optical state ever made 
has a variance of only 7 % of the vacuum variance [2^ . 

Producing squeezed states of itinerant modes allows 
the generation of distributable entanglement by sending 
two copies of a squeezed vacuum state through the two 
input ports of a balanced beam splitter. The coherent 
information [27[ is one useful way to characterize the en- 
tanglement between the two output modes. The asymp- 
totic number of maximally entangled qubit pairs (e-bits) 
that can be distilled per copy of the noisy entangled state, 
by using local operations and one-way classical commu- 
nication, is at least as large as the coherent information 
[2^ . Given two copies of p, one could make two entangled 
modes with 2.5!to4 e-bits of coherent information. 

In conclusion, we have reconstructed the Wigner func- 
tion of an itinerant squeezed microwave field generated 
at the output of a Joscphson Parametric Amplifier. Us- 
ing a second JPA as a preamplifier has increased the 
quantum efficiency of the microwave homodync detection 
from approximately 2 % to 36 %. The level of squeez- 
ing is primarily limited by noise added to the squeezed 
state by the JPA. Improving the performance of the JPAs 
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SUPPLEMENTARY MATERIALS FOR "QUANTUM STATE TOMOGRAPHY OF AN ITINERANT 

SQUEEZED MICROWAVE FIELD" 

DATA ACQUISITION AND CALIBRATION 

Determining the amplifier added noise and loss requires several calibration steps that permit us to isolate the effect 
of a specific loss or added-noise contribution to the overall efficiency of the homodyne measurement. The crucial aspect 
that makes this calibration possible is that the JPA cavities have widely-tunable resonance frequencies, adjusted by 
imposing a magnetic flux Isoj . Far from resonance the JPA cavities behave as open circuits. They are simply 
mirrors that reflect the microwave field without otherwise transforming it; therefore, either the SQ or AMP or both 
stages can effectively be bypassed. 

We begin with both JPA stages bypassed, so that they have Gs — Ga = 1. If the switch were lossless, when 
it is connected to the cold load, the noise power exiting the HEMT amplifier would be 5 = Gh{Ah + Sj), where 
Sf = (1/2) + Uf = (1/2) + [exp{huj/kBTf) — and T/ is the refrigerator's temperature. Notice that the result 
doesn't depend on the transmissivities a, /?, or ^ because these are at the same temperature as the cold load, 
consequently each loss component emits as much power as it absorbs. However, with the switch connected to the 
hot load, the expression for the total power at the output becomes S = Gh{Ah + {£,a(3)Sh + (1 — ^af3)Sf), with 
Sh ~ (1/2) + nil ~ (1/2) + [exp{Uuj/kBTh) — 1]^^ and Th = 4.1 K. In both cases, we expect and observe that S 
depends linearly on Sf with an offset. By fitting these linear dependencies we can extract Gh, Ah, and the product 

We cannot assume that the switch is lossless. Because its loss sits at 4.1 K, it will always emit noise power 
Sh{y — A) + S'i„A, where Sin is the incident noise and A is the switch transmissivity. So, even when n/ <C 1/2, the 
state presented at the SQ stage will have average thermal occupancy n = (1 — X)£^Tih. We write the noise power at 
the output as function oi Sf, for switches in both positions as 

Sic = GhAh + ShGnil ~ X)^al3 + SfiGnX^aP + Gh{1~ ^a/S)] - + mi^Sf (2) 



Sih = GhAh + ShGni^ap) + 5/[Gh(1 - ^ap)] = bih + mihSf, (3) 

where the subscript Ic (l/i) corresponds to the switch connected to the cold (hot) load. Fitting our noise data to 
the right hand side of Eq. [2] and [3l we can obtain the four parameters bih,bic,mih and mic- However as these 
parameters are not independent, Sh = {bih — bic)/{mic — mih), we cannot extract the switch loss independently. We 
can nevertheless bound this unknown loss by taking a worst case estimate as the manufacturers minimum specified 
transmission (at room temperature) A = 0.83 and assuming it is less lossy at 4.1 K. We moreover confirmed that 
at room temperature the frequency dependent loss of the switch is within the manufacturer's specification. Then 
by using 1 < A < 0.83, we can bound the desired parameters using Eq. [2] and [3l with the expressions (^a/?)^^ = 
1 + mihShX/ibih - bic), and Gh = mih/il - ^ap), and Ah = {blc/Gn) - (1 - A)S',,(^a/3). 

We then perform the same analysis, finding the linear dependence of the output noise on Sf and on the switch 
position, with the AMP ON and SQ OFF. From these fits and knowledge of Ah and Gh we find ^a, Aa, andG^/S. 
Finally, we operate the experiment with AMP OFF and SQ ON. A third time we fit the linear dependence of S on Sf 
with the switch in both positions, determining ^, a and /? separately (Fig. [S]d). We evaluate the expressions for a, /3, 
Aa, Ah, Gh and n at the bounds on A, finding the range of values in the main text. We also find ^ = —9.9 ± 1 dB, 
of which 6 dB arises from an attenuator that has been placed at the input of the SQ stage. 

To acquire these calibration data sets, we regulate the refrigerator's temperature at 10 values between base temper- 
ature (T < 50 mK) and 800 niK, which requires about 7 hours to complete. For each temperature point we measure 
the noise at the output under all six conditions, 2 switch positions, and 3 amplifier configurations (AMP OFF SQ 
OFF, AMP ON SQ OFF, and AMP OFF SQ ON). We inject a tone detuned from the AMP pump by 20 kHz. By 
dividing the noise power at the output of the chain by the power in this tone, we become insensitive to any variation 
in Gh over the time needed to acquire the data. At the end of the calibration, we immediately operate the experiment 
with SQ ON and AMP ON, to acquire the data in the paper. In addition, we use the tone to ensure that we do not 
saturate the amplifier chain. 

Data is acquired by digitizing the output (IF port) of the mixer at rate of 10*" samples per second. We filter the IF 
port with a 5 MHz anti-aliasing low-pass filter. The digitized data is digitally filtered with a 3rd-order Butterworth 
high-pass filter with a 500 kHz corner (3 dB) frequency. The noise density S is the average noise density in the 
frequency range between 500 kHz and 5 MHz. 
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FIG. 5: The noise density S in arbitrary units at the output of the measurement versus refrigerator temperature T/. a.) Data 
acquired with the AMP and SQ OFF and the switch connected to the cold load (circles) and hot load (squares). The lines are 
linear fits to S versus 5*/ for the case of the switch connected to the cold load (solid) and hot load (dashed), b.) Data acquired 
with the AMP ON and SQ OFF (blue) and AMP OFF and SQ ON (red), with the switch connected to the cold load (circles) 
with a linear fit (solid) and hot load (squares) with a linear fit (dashed). The arbitrary y-scale is consistent between the six 
plots. The linear fits do not appear as lines because we plot S versus Tf rather than Sf. 



TABLE I: Inferred properties of the squeezed state, upon our three analysis assumptions. 

Pessimistic Best guess Optimistic 

Fidelity 0.66 ± 0.02 0.807 ± 0.016 0.960 ± 0.005 
Min. var. of comparison pure state |^/')": 

0.065 ± 0.009 0.060 ± 0.003 0.0493 ± 0.0006 

p's purity 0.62 ± 0.02 0.74 ± 0.02 0.96 ± 0.01 

p'ssq. var.' 0.918 ± 0.002 0.484 ± 0.013 0.304 ± 0.008 

p's anti-sq. var.*" 25.54 ± 0.07 20.17 ± 0.06 19.18 ± 0.05 

Coherent info. 2.19 ± 0.08 2.46 ± 0.09 3.42 ± 0.05 

Linear sq. var.** 0.40 ± 0.02 0.12 ±0.02 -0.18 ±0.02 

"Ratio of the variance of the squeezed quadrature of the pure squeezed vacuum state with highest fidelity to the variance of the vacuum. 
'Ratio of variance of most likely state p's squeezed or anti-squeezed quadrature to the variance of the vacuum. 
'^Coherent information (in e-bits) that could be produced with two copies of the squeezed state and a beam splitter. 
''Direct linear inference of the squeezed state's minimum variance, relative to vacuum variance. 



MAXIMUM LIKELIHOOD ANALYSIS OF THE SQUEEZED STATE 

Tabic H] shows the statistical errors in our estimates of inferred parameters characterizing the squeezed state, for 
the three analysis cases, based upon our systematic calibration uncertainties. The first line presents the fidelity 
F = {ip\p\ip), where p is the maximum likelihood reconstructed density matrix of the field exiting the SQ and is 
the pure vacuum squeezed state that maximizes the fidelity. The second line gives the ratio of the minimum variance 
of to the variance of vacuum. The third line gives the purity Tr(p^) of p. The fourth and fifth lines give the 
ratios of the squeezed and anti-squeezed variances of the reconstructed state to the variance of vacuum. The sixth 
line presents the coherent information that could be obtained by combining on a beam splitter two copies p. The last 
line gives our estimate the the experimental states' minimum variance based on direct linear inference. 
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We have stated three variances that characterize the state created in this experiment: the hnear estimate of the 
experimental state's minimum variance (12 %), the most hkely state p's minimum variance (48 %), and the minimum 
variance of the pure squeezed vacuum state \ip) that maximizes the fidehty with p (6.0 %). Here we give more 
discussion of these variances. 

The quadrature measurements we observe are the hnear combination of the quantum state created by the squeezer 
and vacuum fluctuations: 



Xe = y/rjx9 + V(l - 'n)ye, 

where xe is the quadrature of the squeezed state, and ye is the quadrature of the vacuum state. Solving for xe gives 

xe = ^ (^Xe - V(l-??)ye) • 

Therefore the inferred variance of the squeezed state's quadrature Axg is 

Axj^^[AX^-{l-v)Ayj)]. 

The vacuum variance Ayg = 1/2, and we can easily calculate an unbiased estimate of AXg for every phase 9. This 
gives us an unbiased estimate of Axg that does not depend on the details (for example, Gaussianity) of the quantum 
state. We calculate Axg using 20,000 quadrature measurements at each of 100 evenly spaced 6 and calculate the 
minimum value Aa;gQ]y[jj^, in Table HI The statistical uncertainties show one standard deviation in the estimate of 
Axgq MiN- ^'-'^ "optimistic case" we calculate a negative variance, which is clearly unphysical. This is a sign of 
inconsistency in the "optimistic" calibration parameters. Because the "optimistic" estimate for the squeezed state is 
computed using the lower bounds on rj and n, this negative variance is evidence that the detector's true rj and / or 
effective gain (AXgg off/^^s^q,Off) must be larger than the lower bounds set by calibration. 

The minimum variance of p is significantly higher than this linear estimate. This is caused by bias in the maximum 
likelihood method. Quantum state estimation by maximum likelihood is biased toward more mixed states, and the 
amount of bias increases with increasing purity of the state from which the measurements arc drawn [31|. Based 
on numerical experiments, the bias in our estimates of the fidelity should be well below the uncertainty level set by 
systematic effects. However, the bias in our estimates of the minimum variance of the inferred state could be larger. 
To attempt to quantify this effect, we simulated measuring and performing maximum likelihood tomography on a 
Gaussian state. This Gaussian state is chosen to have minimum and maximum variances equal to those calculated 
by the linear method described above for the "best-guess" case. By computer we simulate 10,000 quadrature mea- 
surements (the same number we used for ML analysis of the true experiment) from this Gaussian state and perform 
maximum likelihood tomography on those measurements. The inferred state has minimum variance 40 %. Therefore 
it is possible that the experimental state has smaller minimum variance than the most likely state inferred from only 
10,000 measurements. Because we have some independent evidence for non-Gaussian effects in the experiment, we 
cannot quantify this size of this bias using this Gaussian simulation. Other numerical simulations have confirmed that 
this bias decreases as the number of measurements analyzed increases and that this bias is not caused by truncation 
of the Hilbert space at 30 photons. 

The apparent discrepancy between the 6 % for the variance of \4>) and the 48 % for the variance of p also deserves 
some comments. It is important to note that one would not expect the minimum variance of a mixed state to equal 
the minimum variance of its highest fidelity pure state. The fidelity between a mixed Gaussian state (centered at the 
origin of phase space) whose minimum and maximum variances are Vx and Vp and a pure squeezed vacuum state with 
minimum variance Vs is given by 



F, 



Gauss 



The highest fidelity pure state has minimum variance Vs = \f^^ ^^"^ ^^"^ fidelity between these two states is 

2 

-^Gauss.max — "j ; , 

1 + 2^VxVp 

Consider the state cr to be a Gaussian state with minimum variance of 48 % and maximum variance 2017 %. (cr 
has variances equal to those of our state p, but unlike p, a is guaranteed to be Gaussian.) Then let 1-0) be the pure 
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squeezed vacuum state that has maximum fidehty with <t. i^Causs.max = = 0.49, and the minimum variance of 

IV') is 7.7 %. The difference between the minimum variances of p and \'4)) is to be expected. However, the maximum 
fidehty of p is significantly larger than we would expect if it was perfectly Gaussian. This non-Gaussianity could be 
caused by bias in the maximum likelihood inference and / or genuine non-Gaussian effects in the experiment. 

Tomographic reconstruction of a quantum state requires that the experimental device always creates the same 
(potentially mixed) quantum state, that the measurements are well described by inefficient quadrature measurements, 
and that the calibration of those measurements is consistent. In this experiment we have observed some evidence that 
at least one of these assumptions is violated. The likelihood of the maximum likelihood state is significantly lower 
than one should expect from simulated measurements on that state. That is, if the tomographic assumptions above 
were true, we expect to find a significantly higher value for the maximum likelihood. We believe this effect could be 
caused by an interaction between the state preparation and measurement stages of the experiment, such as a phase 
dependent efficiency of the measurement JPA, and/or non linear processes in the measurement. 



